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The "finite-size" effects in isobaric yield ratio (lYR), which is shown in the standard grand- 
canonical and canonical statistical ensembles (SGC/CSE) method, is found to prevent obtaining 
the actual values of physical parameters. Due to there is no inputs of neutron and proton density 
distribution in SGC/CSE, the conclusion of SGC/CSE is doubtable. To investigate wether lYR has 
the "finite-size" effects, lYR of mirror nuclei are reexamined using a modified statistical abrasion- 
ablation model. lYR of mirror nuclei in the UOA MeV ^'^^'^^Ca + ^^C, and ^^Zn/^^Zr/^^^Sn + ^^C 
reactions are calculated using a modified statistical abrasion-ablation model. lYR of mirror nucleus 
in reactions of ^*~^^Ca increases when the projectile becomes more neutron-rich, but saturates in 
the reactions of very neutron-rich nucleus. lYR of mirror nuclei in reactions of ^^Zn/^^Zr/^^^Sn 
that have same n/p ratio show no "finite-size" effects. lYR of mirror nuclei is found to depends on 
the isospin of the projectile rather than the size of reaction system. 

PACS numbers: 25.70.Pq, 21.65.Cd, 25.70.Mn, 24.60.-k 



I. INTRODUCTION 

The yield of fragment is often constructed to extract 
information of nuclear matter and temperature in heavy- 
ion collisions (HIC). The double yield ratio of light frag- 
ments are constructed as thermometer [iHl] . The scaling 
of the ratios between isotopic yield in similar reactions, 
i.e., the isoscaling, is used to study the symmetry energy 
of hot emitting source in framework of various theoretical 
models IMS. 

Using a modified Fisher model (MFM) that based on 
free energy [3: [3 > the isobaric yield ratio (lYR) is used 
to extract the symmetry energy of fragment in HIC near 
the Fermi energy The correlations between lYR 

and energy terms which contribute to the free energy of 
fragment are constructed to extract physical parameters 
associated with mass of nucleus [2l|. In lYR methods, 
the energy terms that only depend on mass cancel out, 
which make it convenient to investigate the symmetry- 
energy term in mass formula. 

Recently, it is discussed that the lYR method, which 
utilizes only one reaction system, does not provide can- 
celation or minimization of the effects associated with 
mass and charge constraints due to the "finite-size" ef- 
fects in the standard grand-canonical and canonical sta- 
tistical ensembles (SGC/CSE) 

The linear behavior of lYR for the mirror nuclei, which 
is discussed in Ref. (20l - [2^ . is written as, 

IYR{m) =ln[i?(l,-l,A)] = Y{A,l)/Y{A,~l) 

= [Mn - + 2ac(^ - l)/A^/^]/T, (1) 
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where ac is the Coulomb-energy coefficient, ^„ and 
are the chemical potential of neutron and proton, respec- 
tively. 

The "finite-size" effects of lYR(m) in SGC/CSE, which 
depend on the volume of reaction system, is discussed to 
prevent one from obtaining precise information on the 
nuclear properties using the isobaric ratios method. It's 
also pointed out that the "finite size" effects are negligi- 
ble only for system sizes much larger than those actually 
formed in the experiments [l^. Actually, though the fi- 
nite effects are believed to be minimized from the n/p 
ratio fixed source in isoscaling, the finite size effects are 
still obvious in isoscaling in some models [Ti- 
lt's well known the yield of fragment has strong depen- 
dence on the isospin of the projectile, and the neutron 
and proton density distributions inside the projectile. 
Actually, the neutron and proton density distributions 
arc important inputs in nuclear reaction models, while 
there is no such information in SGC/CSE [22,[23]- The 
temperature has a wide distribution due to the strength 
of collisions, while it is set to a certain value in SGC / CSE. 
The yiel ds of fragments in SGC/CSE are not shown in 
Ref. [22| ■ Considering the density distributions of nuclei 
involved, SGC/CSE is not a good model to evaluate the 
yield of fragment, especially in reactions of neutron-rich 
nucleus. Thus the deduction based on yields of fragments 
in SGC / CSE may be doubtable. The statistical abrasion- 
ablation (SAA) model can well reproduce the yield of 
fragments in HIC induced by neutron-rich nucleus at in- 
termediate energy and relativistic energy 0, H^H^ , and 
it is used to study the isospin effect in HIC. Consider- 
ing the dependence of yields and the resultant lYR(m) 
on the isospin of projectile, the "finite-size" effects in 
lYR(m) shown in SGC/CSE (which depend on the vol- 
ume of reaction system), can also be explained by the 
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isospin effect in isotopic yield. 

In this article, using the SAA model, the lAOA 

MeV ^^'^^^Ca which have different isospins, and 
■^^Zn/s^Zr/^^ogj^ l-l^j^l- i^j^g gg^j^g j^/p j.g|-jQ 48cj^ (which 

have been calculated in Ref. [22|) are calculated to show 
that lYR(ni) depend more on the isospin but not the size 
of the projectile, i.e., the "finite-size" effects of lYR(m) 
in SGC/CSE should be explained as the isospin effect. 



II. MODEL DESCRIPTION 

The SAA model was developed by Brohm and Schmidt 
to predict isotopic yield in the relativistic energy HIC [l^, 

fand modified to describe the intermediate energy HIC 
m. Hi, It can weU reproduce the yields of 

both the small and large mass fragments in the lAOA 
MeV 40.48Cjj_^9Bg g^j^jj 58,64]s^i_^9Bg projectile fragmen- 
tation j25l [3l| . For simplification, only some important 
formula are listed here since the model is well described 
inRefs. [23, [11,123,133, p. 

The SAA model takes independent nucleus nucleus col- 
lisions for participants in an overlap zone of the two col- 
liding nuclei and determines the distributions of abraded 
neutrons and protons. The colliding nuclei are described 
to be composed of parallel tubes orienting along the beam 
direction. Neglecting the transverse motion, the collision 
is described by independent interactions of tube pairs. 
Assuming a binomial distribution for the absorbed pro- 
jectile neutrons and protons in the interaction of a spe- 
cific pair of tubes, the distributions of the total abraded 
neutrons and protons can be determined. For an in- 
finitesimal tube in the projectile, the transmission prob- 
abilities for neutrons (protons) at a given impact param- 
eter b are calculated by 

tk{s- h) = exp{-[Z?^(s - 6)a„fc + - b)a.pk\] (2) 

where D"^ is the nuclear-density distribution of the tar- 
get integrated along the beam direction and normalized 
by ^(PsDl = N'^ and /d^si^J = Z^. With N'^ and 
.Z"^ refer to the neutron and proton number in the tar- 
get respectively, the vectors s and b are defined in the 
plane perpendicular to beam, and ayk is the free nucleon- 
nucleon reaction cross section. The average absorbed 
mass in the limit to infinitesimal tubes at a given b is 

< AA(&) j (fsDl{s)[l-t^(s-b)] 

+ j SsDP{s)[l-^{s-b)] (3) 

The proton and neutron density distributions are as- 
sumed to be the Fermi-type, 



i + cxp(f7ri)^ 



n,_p 



(4) 



where pf is the normalization constant, ti is the diffuse- 
ness parameter, and Ci is half the density radius of the 



neutron or proton density distribution. The parameters 
in Eq. ^ can be different for proton and neutron in 
asymmetric nucleus. 

The cross section for a specific isotope can be calcu- 
lated from 



ct(AA^,AZ)= J (fbP{AN,b)P{AZ,b), 



(5) 



where P{AN, b) and P{AZ, b) are the probability distri- 
butions for the abraded neutrons and protons at a given 
impact parameter b, respectively. 

The second stage of the reaction is characterized by 
the evaporation of particles, of which in the SAA model 
is described by the conventional statistical model un- 
der the assumption of thermal equilibrium in the excited 
prefragment. The excitation energy of projectile spec- 
tator is estimated by a simple relation of E* = 13.3 < 
AA{b) >MeV, where 13.3 is the mean excitation en ergy 
due to an abraded nucleon from the initial projectile [27|. 
After de-excitation, isotopic yield comparable to experi- 
mental result is obtained. 
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FIG. 1. 



(Color online) Tfie isobaric yield ratio for mirror 
nuclei lYR(m) in the 140 A MeV ^^^^^Ca + ^^C reactions 
calculated by the SAA model. 



III. RESULTS AND DISCUSSIONS 

The isotopic yields in the UOA MeV ^^^^'^Cn + ^^C, 
and ''^Zn/^^Zr/^^^Sn+^^C reactions are calculated us- 
ing the SAA model. lYR(m) of the calculated 140 A 
MeV 38~52(2;a^ _|„ 12q reactions are plotted in Fig. [H It 
is found that lYR(m) in the '^^'^^^Ca reactions increase 
when the projectile nucleus becomes more neutron-rich. 
Defining x = 2{Z — 1)/A^^^, when x < 6, the difference 
among lYR(m) of these reactions are quite small. Fairly 
well linear correlations between lYR(m) and x are found 
in the ^^'^"^Ca reactions. For the neutron-deficient '^'^Ca 
and neutron-proton symmetric ^'^Ca projectiles, lYR(m) 
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firstly increase linearly with x when a; < 7, but de- 
crease when X > 7 . For the neutron-rich "'^'^^^Ca pro- 
jectiles, lYR(m) increase linearly when x are not large, 
but lYR(m) increase faster when x of fragments are lager 
than 0.85x of projectile. For these fragments, only few 
nucleons are removed from the projectile, which are very 
typical projectile-like fragments. It has been illustrated 
that these fragments arc mostly produced in the periph- 
eral collisions with large impact parameters in the SAA 
model j33j . The mirror nuclei with small x and large mass 
X have different production mechanisms. Considering the 
slopes of the linear correlations [for the linearly increas- 
ing part in lYR(m)], the slope increases as neutrons of 
the projectile increases. For ^''Ca and ^^Ca, lYR(m) al- 
most overlap except few fragments with large cc, which 
shows the saturation of lYR(m) in reactions induced by 
the very neutron-rich nuclei. 



10 



i 5 

£ 4 

> 3 



* Ca * Ca exp 
. '°Ca e '"Caexp 



°Zr 




12 



14 



x=2(Z-1)/A" 



FIG. 2. (Color online) The isobaric yield ratio for mirror nu- 
clei. The full and open circles (stars) are for the SAA and 
experimental 140^4 MeV ^''Ca(**Ca) induced reactions, re- 
spectively; the triangles, hexagon, and diamonds represent 
lYR(m) for the 140A MeV ■^^Zn/^'^Zr/^^^Sn + ^^C reactions. 
The lines are the linear fitting results by Eq. 



In Fig. [21 the SAA calculated and experimental 
lYR(m) in the 140 A MeV 40,48^3^ induced reactions 
are plotted as the solid and open symbols, respectively. 
lYR(m) in the 140 A MeV ^^Zn/^^Zr/i^OSn H- ^^G re- 
actions are also plotted. Though the experimental odd- 
even staggering in lYR(m) can not be reproduced, the 
SAA results fit the ^"^Ca experimental data rather well, 
including the decreases of lYR(m) when a; > 9. For 
^^Ca, the SAA lYR(m) are smaller than the experimen- 
tal ones. The linear fitting of the SAA and experimen- 
tal lYR(m) read y = (0.62 ± 0.02)a; - (1.24 ± 0.16) and 
y = (0.75 ± 0.06)x - (0.92 ± 0.50), respectively. The 
large difference among lYR(m) in reactions induced by 
projectiles with same n/p ratio in SGC/CSE disappears 
in the SAA results. lYR(m) for ^^Zn projectile exhibits 
quite good linear correlation when x < 11, which shows 



very little difference to the lYR(m) for ''^Ca. The smah 
mass mirror nuclei do not survive the deexcitation pro- 
cess in SAA calculation for ^^Zr and ^^OSn, but lYR(m) 
for ^^Zn, ^^Zr, and ^^°Sn show very little difference when 
X < 12. Considering the saturation of lYR(m) in large 
isospin nucleus induced reactions, the overlap of lYR(m) 
of ^^Zn, ^^Zr, and ^^°Sn can be explained. 
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FIG. 3. (Color online) a^/T and fi/T fitted from lYR(m) by 
Eq. (HJ. The labels ""Ca*, "^Ca*, and stars represent the 
values for the experimental results. 

In Fig. m the fitted ac/T and A^i from lYR(m) in 
the ^S'-'^sca and ^^Zn/^^Zr/nzgn reactions are plotted. 
The fittings are limited to the linear increasing part of 
lYR(m). flc/T (A/i/r) of lYR(m) increases (decreases) 
when the projectile becomes more neutron-rich, but ac/T 
(A/x/r) are similar due to the similar lYR(m) in reac- 
tions of the neutron-rich projectiles. For projectile with 
the same n/p, Qc/T (A/u/T) has very little difference and 
no finite size effects are shown. It's easy to conclude that 
lYR(m) depends more on the isospin of the projectile 
rather than the size of system, and the isospin depen- 
dence of ttc/T (A/i/r) is weakened and even disappears 
in reactions of neutron-rich nucleus. 

According to Eq. (O, the yield of a prcfragmcnt is 
mainly decided by the nucleon-nucleon cross sections and 
the proton and neutron density distributions. Differ- 
ently to the strictly linear increase of lYR(m) with x in 
SGC/CSE, lYR(m) of fragments with small x has very 
little difference but show accelerating changes in frag- 
ments with large x in SAA. This phenomena can be ex- 
plained as the isospin effect in HIC, i.e., the similarity in 
lYR(m) of small x due to the similarity of proton and 
neutron density distributions in central collisions, while 
the large difference in lYR(m) of large x are correspond- 
ing to the large difference between the proton and neu- 
tron density distributions 0, [Ml, [H, [s^. The obvious 
lYR(m) dependence on the isospin of the projectile also 
exhibits the behavior of the " finite size" effects shown in 
SGC/CSE. But contrary to the the "finite size" effects 
in SGC/CSE, the dependence of lYR(m) on the system 
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volume of projectile with same n/p is found to be very 
small. 

Using one reaction, the extracted /ip)/T and Oc/T 
using the lYR method may depend on the n/p ratio of the 
projectile. In Huang's workj2(|. the Cc/T is scaled to the 
reaction system parameter Z/A(= {Zp + Zt)/{Ap + At)). 
If using the methods in Ref. |2C| to extract the cLsym/T 
of fragments, an increase of 0.15 of ac/T may introduce 
a difference of 0.6^1.5 when x increase from 4 to 10, 
i.e., for fragment with small a;, the difference is rather 
small. Based on these results, it is proposed the isobaric 
yield ratio has very little "finite-size" effects depending 
on the size of reaction system, and can be used to extract 
the coefficients of energy terms contributing to the free 
energy. 

To summarize, the "finite size" effects in lYR(m) in 
SGC/CSE are revisited using a modified SAA model. 
The SAA results well reproduce the lYR(m) in the mea- 
sured lAQA MeV '*°Ca-|-^Be reactions, while the SAA cal- 
culated lYR(m) is smaller than the experimental results 
in the 140 A MeV "^Ca-H^Be reaction. lYR(m) is found 
to depend on the isospin of projectile nuclei. lYR(m) 



saturates and its dependence on the isospin of projec- 
tile is weakened and even disappears in reaction of very 
neutron- rich nucleus. The "finite-size" effect shown in 
GC/SCE disappears in the SAA resuhs. The lYR(m) 
depends on the isospin of the projectile rather than the 
size of reaction system shown in the results of SGC/CSE. 
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